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This article presents a phenomenological dynamic phase transition theory - modeling and analysis 
- for liquid helium-3. We derived two new models, for liquid helium-3 with or without applied field, 
by introducing three wave functions and using a unified dynamical Ginzburg-Landau model. The 
analysis of these new models leads to predictions of existence of 1) a unstable region, 2) a new 
phase C in a narrow region, and 3) switch points of transition types on the coexistence curves near 
two triple points. It is hoped that these predictions will be useful for designing better physical 
experiments and lead to better understanding of the physical mechanism of superfluidity. 

PACS numbers: 



I. INTRODUCTION 

Superfluidity is a phase of matter in which "unusual" 
effects are observed when liquids, typically of helium-4 or 
hclium-3, overcome friction by surface interaction when 
at a stage, known as the "lambda point" for helium-4, at 
which the liquid's viscosity becomes zero. Experiments 
have indicated that helium atoms have two stable iso- 
topes 4 He and 3 He. 3 He contains two electrons, two pro- 
tons and one neutron. Hence it has a fractional spin and 
obey the Fermi-Dirac statistics. The liquid 3 He has two 
types of superfluid phases: phase A and phase B. In par- 
ticular, if we apply a magnetic field on the liquid 3 He, 
then there will be a third superfluid phase, called the A\ 
phase. 

The main objectives of this article are 1) to establish 
some dynamical Ginzburg-Landau models for 3 He with 
or without applied magnetic fields, and 2) to study super- 
fluid dynamic transitions and their physical significance. 

Consider the case without applied magnetic field. The 
modeling is based on two main ingredients as follows. 

First, instead of using a single wave function tp, we use 
three wave functions to represent Anderson-Brinkman- 
Morel (ABM) state and the Balian-Werthamer (BM) 
state. More precisely, we introduce three complex val- 
ued functions ipQ,ip\, ip2 to characterize the superfluidity 
of 3 He , with tp for the state | |T>, ipi f° r the state | H>, 
and vj 2 for the state | ti> +| IT>- Then we are able to 
formulate a Ginzburg-Landau (GL) energy in terms of 
these three wave functions and the density function p n 
of the normal fluid state. 

Second, we use a unified time-dependent Ginzburg- 
Landau model for equilibrium phase transitions, devel- 
oped recently by the authors [4j [5], to derive a general 
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FIG. 1: The coexistence curve of 3 He in the case without an 
applied magnetic field. 

time-dependent GL model for 3 He. 

In a nutshell, the model is obtained by a careful ex- 
amination of the classical phase transition diagrams and 
by using both mathematical and physical insights offered 
by a recently developed dynamical transition theory as 
briefly recalled in the appendix. The model of the case 
with applied magnetic field can be derived in the same 
fashion. 

With the models in our disposal, we can study the 
dynamic phase transitions of liquid 3 He, and derive some 
physical predictions. 

To be precise, we first recall the classicalphase tran- 
sition diagrams of 3 He, as shown in Figure [I] see among 
many others Ginzburg [T], Reichl [7] and Onuki [5J. 
From this diagram, we see that there are two coexis- 
tence curves, with one separating the solid state and the 
phase A and phase B superfluid states, and the other 
separating the superfluid states and the normal liquid 
state. The transition crossing the coexistence curve be- 
tween the solid and superfluid states is first order (Type- 
II in the sense of dynamic classification scheme given 
in the appendix), and the transition crossing the coex- 
istence curve between the superfluid and normal liquid 
states are second order (Type-I in the sense of dynamic 
classification scheme) . In addition, the second transition 
between phase A and phase B superfluid states is first 
order (Type-II). 

We would like to mention that the Ginburg-Landau 
theory with only one wave function can not describe this 
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phase transition diagram, and this is one of the mains 
reasons that we need a new Ginzburg-Landau model as 
discussed above to study the phase transition dynamics 
for 3 Hc. 

The models established are analyzed using a recently 
developed dynamical transition theory, leading to some 
interesting physical predictions. Here we address briefly 
the new results derived. For the case without the applied 
magnetic field, the main results obtained are synthesized 
in a theoretical PT-phase diagram of 3 He given by Fig- 
ure n 
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FIG. 2: Derived theoretical PT-phase diagram of 3 He: The 
region H = Hi U H2 is the unstable domain, where the solid 
state and the superfluid state appear randomly depending on 
fluctuations. The curve bed is the first critical curve where 
phase transition between normal fluid and superfluid states 
occur. 

One prediction from our results is the existence of a 
unstable region H — Hi U H2, in which the solid state 
and the superfluid states A and B appear randomly de- 
pending on fluctuations. In particular, in Hi, phase B 
superfluid state and the solid may appear, and in H2 the 
phase A superfluid state and the solid sate may appear. 

Another prediction is the possible existence of phase 
C superfluid state, which is characterized by the wave 
function ip 2 , representing | fj,> +| J,f>. However, phase 
C region is very narrow, which may be the reason why it 
is hard to be observed in experiments. 

Also, the results predict that near the two triple points 
b and c, there is a possibility of the existence of two switch 
points, where the transition on the corresponding coex- 
istence curve switches types at each switch point. The 
existence of such switch points depends on the physical 
parameters. 

One important new ingredient for the analysis is a dy- 
namic transition theory developed recently by the au- 
thors [2J [3] . With this theory, we derive a new dynamic 
phase transition classification scheme, which classifies 
phase transitions into three categories: Type-I, Type-II 
and Type-Ill, corresponding respectively to the contin- 
uous, the jump and mixed transitions in the dynamic 
transition theory. 

The case with applied field can be addressed in the 
similar fashion; see Section V. 

This article is organized as follows. Section II intro- 
duces a new dynamic model for liquid helium-m without 



applied magnetic field. The phase transition dynamics is 
given in Section III, and the type of transitions in dif- 
ferent regimes are given in Section IV. Section V deals 
with the case with applied magnetic field, and Section VI 
gives a summary on the physical predications of the new 
models. 



II. 



DYNAMIC MODEL FOR LIQUID 3 HE WITH 
ZERO APPLIED MAGNETIC FIELD 



The superfluidity of liquid 3 He was found in 1971 by 
D. M. Lee, D. D. Osheroff, and R. C. Richardson, and 
its transition temperature is T ~ 10 _3 -ftT under p = 1 
atm (10 5 Pa). There are two superfluid phases A and B 
if there is no applied magnetic field. Figure [T] provides 
the phase diagram of the liquid 3 He in the PT-planc. 

Because the atoms 3 He are fermions, to form the su- 
perfluid phase they must be paired to become bosons. 
When no magnetic field is applied, there are two pairing 
states: the Anderson-Brinkman-Morel (ABM) state and 
the Balian-Werthamer (BW) state, given respectively by 

V^l* >= (a + ib)\ |T> +(a-ib)\ U>, 

%/2> >= (a + ib)\ TT> +(a - ib)\ U> +c[\ U> +1 IT>]. 

The ABM state corresponds to the superfluid phase A, 
and the BW state to the phase B. Also, we call state 
I T4> +1 4T> as phase C, which appears in the theory 
developed in this article. 

We introduce three complex valued functions 
ipo,ipi,ip2 t° characterize the superfluidity of 3 He , 
in which ipg to the state | t|>,Vi to the state | H>, 
and 1P2 to the state | Ti> +| jt>. 

Let p n be the normal fluid density, p a , Pb and p c rep- 
resent the densities of superfluid phases A, B and C re- 
spectively. Then we have 

Pa = Toluol 2 + Ti\ipi\ 2 (t o >0,ti>0), 
Pb = T 2 \tpo\ 2 + r 3 \tfji\ 2 + T 4 \ip 2 \ 2 (r 2 ,r 3 > 0,r 4 > 0), 

Pc = T 5 \^ 2 \ 2 (r 5 >0). 



The total density of 3 He is given by 



Pn 




in 


the normal state, 


Pn ~ 


V Pa 


in 


the phase A state 


Pn ~ 


V Pb 


in 


the phase B state 


Pn ~ 




in 


the phase C state 



Physically, the states V'OjV'i an d ^2 are independent, 
and consequently there are no coupling terms \V(ipi + 
-0j)| 2 and \ipi + ?pj\ 2k (i ^ j) in the free energy density. 
Since in the case without applied magnetic field ipo an d 
1^1 are equal, their coefficients in the free energy should be 
the same. Thus, we have the following Ginzburg-Landau 
free energy for 3 He with H = 0: 
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G(i> ,1pi,1p2,p n ) — 

1 -in 



-|V^o| 2 + ail^ol 2 + a 2 p„|^ | 2 + y IV>o| 4 



-|V^i| 2 + «i|Vi|' 



i , ,9 a 3 i , ,4 k 2 h 2 



2 m 

+ /WlH^I 2 + ^Pnl^P + fc 3 |Vp„| 2 + pi\ Pn 



2 . 2 p 2 3 P3 3 / . MO 2 

+ -g-pn + yPri - P [Pn + y P„ 



da:, (1) 



where the coefficients depends on T and p, and for 1 < 
* < 3, j = 2,3,4, 



> 0, & > 0, 
a 2 > 0, a 3 > 0, p 3 > 0, p 2 < 0. 



(2) 



For ai,P\ and pi, there are regions Ai,Bi 7 Ci (i = 1,2) 



in the PT-plane Rl such that A 1 + A 2 = Bi + B 2 



L+, and 



ai = ai(T,p) 
/3i = Pi {T, p) 

m = pi(T,p) 



>o if(T,p)eA!, 

<0 if(T,p)eA 2 , 

>0 if(T,p)eBi, 

<0 if(T,p)eB 2 , 

>0 if^^ed, 

<0 if (T,p)eC 2 . 



(3) 
(4) 
(5) 



r 



It is known that for 3 He , pi — pi(T, p) is not mono- 
tone on T. In fact, at T m = 0.3l8K,p m = 29.31 x 10 5 Pa, 
we have 



/rri \ n dfil(T m ,P m ) 

Pi(T m ,p m )=0, gj, =°) (6) 



where to = (T m ,p m ) is as shown in Figure [l] Near the 
point to the famous Pomeranchuk effect takes place, i.e., 
when pressure increases, the liquid 3 He will absorb heat 
to undergo a transition to solid state. 



By the normalized model ( A6 1 , we infer from (JT|) the 
following time-dependent GL equations for the superflu- 
idity of liquid 3 He : 



dipo 
dt 


hh 2 

Aip 

TO 


- Ot!lp 


- a 2 p n tpo - 


Ps\H 2 ^o- 


o^ol^o, 


dipt 
dt 


TO 


- Ctllpl 


- a 2 Pn^i - 


ftWVi - 




dip 2 
dt 


k 2 h 2 

Aip 2 

TO 


-Pirh 


-PslM^ 


-PM 2 ^2 


- PiPn^2 - P2\lp2\ 2 1p2 



dp n 
dt 



k 3 Ap n - (pi - p p)Pn - P2P„ - P3P 

The nondimensional form of ([7]) can be written as 



2 .3 a 2 I , ,2 "2 L( . ,2 



m 2 



dt 
(hp_ 
~di 
dj>2 

dt 

dp n 

dt 



Aip + AiV>o - aip„0 o - a2|V'2| 2 V'o - a 3 \ip \ 2 ipo, 
- = Atp! + Ai0i - aip„0i - a 2 |V'2| 2 V'i - aal^if^i, 

KiAip 2 + \ 2 ip 2 - bip n ip 2 - b 2 \^ \ 2 ip 2 - &2|0i| 2 02 - b s \4i 2 \ 2 Tp 2l 
K 2 Ap n + A 3 p„ - ci\i> \ 2 - ci|0i| 2 - c 2 \ip 2 \ 2 + c z p 2 n - c 4 p 3 

I 



(7) 



(8) 



where 



Ai 
A 2 
A 3 



ml 

Why 

ml 2 
Wk x 
ml 2 



(ai + a 2 p°), 

(Pi+hp n ), (9) 
(pi - pop - 2p°p 2 - 3(p°) 2 p 3 ), 



and p° is determined by the state state solution of the 
form (0o,^i,^ 2 ,Pn) = (0,0,0,p0) of @. By (§ the 
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coefficients in (|8| satisfy 

a, ; > 0, b t > 0, Cj > VI < i < 3, 1 < j < 4. 
The boundary conditions are the Neumann conditions 
d 



dn 



an. 



(10) 



When p is a constant, the problem rt8b with ( 10 1 can be 



reduced to the following system of ordinary differential 
equations: 



dpi 
dt 

dpi 
dt 

dp n 
dt 



Xipi - aip n pi - a 2 p2Pi - a 3 pi, 
A2P2 - hp n p2 - hpiP2 - hpl, 



X 3 Pn -Clpx- C2P2 + C 3 pl - C4,pn, 



(11) 



where pi = |-0 O | 2 + |^i| 2 and p 2 = |V>2| 2 - 

III. CRITICAL PARAMETER CURVES AND 
PT-PHASE DIAGRAM 

A. Critical parameter curves 

Critical parameter curves in the PT-plane are given by 



h = {(T, P )eR 2 + \ \i(T,p)=Q}, 



1,2,3, 



where A^ = Xi(T,p) are defined by (pj). 

It is clear that the critical parameter curves U are asso- 
ciated with the PT-phase diagram of 3 He . As in the last 
section if we can determine the critical parameter curves 
h (1 < * < 3), then we obtain the PT-phase diagram. 

Phcnomenologically, according to the experimental 
PT-phase diagram (Figure [TJ, the parameter curves U 
(i = 1,2,3) in the PT-plane should be as schematically 
illustrated in Figure [3]ja)-(c). The combination of the di- 
agrams (a)-(c) in Figure [3] gives Figure [4] in which the 
real line bm stands for the^ coexistence curve of the solid 
and liquid phases, and bed for the coexistence curve of 
superfluid and normal liquid phases. 

Now we rigorously examine phase transitions in differ- 
ent regimes determined by the equations and the critical 
parameter curves. 



B. States in the unstable region 

We consider the dynamical properties of transitions 
for (111 in the unstable region. It is clear that at point 



b = (Tb,p b ), 

Xl(T b ,p b ) = \ s (T b ,p b ) =0, X 2 (T bl p b ) <0, (12) 

and the unstable region Hi near the triple point b is 
defined by 

Hi = {(T,p) e R\ I (X 1 ,X 2 ,X 3 )(T,p) = (+,-,+)}. 




FIG. 3: (a) The curve h (Ai = 0), (b) the curve h (A 2 = 0), 
and (c) the curve Z3 (A3 = 0). 




FIG. 4: The curve abch is h, efed is h, and gfbm is Z3 



To study the structure of flows of (111 for (T, p) e H 
it is necessary to consider the equations ( IT]) at the point 
b = [T bl p b ), and by (12 1 which are given by 



dpi 
dt 

dp n 
dt 



We know that 

c 3 > 



= -aip n pi 



a-iPl, 



-Cipi + c 3 p„ 



C4PI, 



for (T,p)dHi. 



(13) 



(14) 



Equations (13) have the following two steady state solu- 
tions: 

z i = (Pi,Pn) = (0,c 3 /a 4 ), 

Z2 = (pi,Pn) = { —a, -a) , 
\a3 / 



C3 



( _ _ L 4ciaic 4 ^ 

2c 4 1 V a 3 cj 



By direct computation, we can prove that the eigenvalues 
of the Jacobian matrices of ( 13 1 at Z\ and Z 2 are neg- 



ative. Hence, Z\ and Z2 are stable equilibrium points 



of (13). Physically, Z\ stands for solid state, and Z2 



for superfluid state. The topological structure of ([13]) is 
schematically illustrated by Figure ^ a), the two regions 
Pi and R2 divided by curve AO in Figure |5jb) are the 
basins of attraction of Z\ and Z2 respectively. 
We note that in Hi, Xi and A 3 are small, i.e., 

0<Ai(T,j») > A 2 (T,p)«l, for (T,p)eH, 
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and ( 13 ) can be considered as a perturbed system of ( 11 1 



Thus, for (T,p) S Hi the system pd[ ) have four steady 
state solutions Zi — Z(T,p) (1 < i < 4) such that 

lim (Zi,Z 2 ,Z 3 ,Z4)(T,p) = (Zi,Z 2 ,0,0), 

(T,p)— (T c ,po) 

and Zi and i? 2 are stable, representing solid state and 
liquid He-3 state respectively, Z 3 and Z4 are two saddle 
points. The topological structure of ( [13) for (T,p) 6 i?i 
is schematically shown in Figure [5^c), and the basins of 
attraction of Z\ and Z2 are i?i and R 2 as illustrated by 
Figure |5jd) . 




FIG. 5: 

C. First phase transition 

On the coexistence curve be, 

Xi = 0, A 2 < 0,A 3 < 0. 

Hence the first phase transition crossing be is between 
normal fluid state and phase A superfluid state. 
On the coexistence curve cd, 

Ai < 0, A 2 = 0,A 3 < 0. 

In this case, the first phase transition crossing this coex- 
istence curve is between the normal fluid state and the 
phase C superfluid state. 

D. Second phase transitions 

When (T,p) crosses the curve segment bed, (l~l| will 
undergo a second transition. We need to consider two 
cases. 

Second transition crossing f'c 



If (T, p) passes through this curve segment be, then the 
first transition solution is given by 

(pl,P2,Pn) = (p*,0,Pn), 

p* 1 =p* 1 (T,p)>0, 
pi = p* n (T,p) < 0. 

Take a transformation 

Pi = Pi- Pi, P2 = P2, Pn=Pn-pn- 



Then, the system (111 is in the following form (for sim- 
plicity, we drop the primes): 



dpi 
dt 

dpi 
dt 

dfn 
dt 

where 



A1P1 - aip\p n - a 2 p\p 2 - a 2 piP2 - a 3 pl 
A 2 p 2 - bip n p 2 - b 2 pip 2 ~ b z p\, 
- A 3 /o n - cipi - c 2 p 2 + (c 3 - 3c 4 p* )p\ 



Ai = Ai + ai\p* n \ - 2a 3 pl, 
A 2 = A 2 + 61 1/5* I - b 2 p\, 
A 3 = A 3 - 2c 3 \p* n \ - 3c 4 p* 2 . 



(15) 



(16) 



The linear operator in (15) reads 




L 



The three eigenvalues of L are 

1 " 



Pi = A 2 , P ± = 



Ai + A 3 ± J (A 3 - Ai) 2 + Acucipl 



(17) 

It is known that the transition solution (p\ , 0, p* n ) is stable 
near be. Therefore the eigenvalues of L satisfy 

Pi(T,p) < 0, (3±(T,p) < for (T,p) near be. 

However, near fc there is a curve segment f'c such that 

fe= {(T,p) ER 2 + \ pi(T,p) = 0}. 



Thus, system ( 15 ) has a transition on f'c, which is called 



the second transition of ([11]), and f'c is the coexistence 
curve of phases A and B; see Figure [2] 

Second transition crossing eh' 

If (T, p) passes through this curve segment be, then the 
first transition solution is given by 

(pl,P2,Pn) = (0,J?2,»7n), 

m = m( T ,p) > °> 

Vn = r) n (T,p) < 0. 
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Take a transformation 



Pl=P, P2=P2~m, Pn = Pn-Vn- 



Then, the system (111 is in the following form (for sim- 
plicity, we drop the primes): 



dpi _ 
dt 

dpi _ 
dt 

dp n _ 

dt 
where 



Aipi - aipip n - a 2 piP2 - a 3 pj, 

X2P2 - hmPn - hmPi - hp2Pn - b 2 pip2 - b 3 p\, 

A 3 p„ - cipi - c 2 p2 + (C3 - 3c 4 ?7„)p^ - c 4 p^, 

(18) 



Ai = Ai + ai|r? n | - a 2 f? 2 , 
A 2 = A 2 + 6i|?7„| - 2b 2 rj 2 , 
A3 = A 3 - 2c 3 |?7„| - 3c 4 ??£. 



(19) 



The linear operator in (18 1 reads 



t I 



L = 







-0 \ 



-6 2 ??2 A 2 -hr]2 
V -ci -c 2 A 3 / 



The three eigenvalues of L are 



ft = Ai, 0± = 



Ai + A 3 ± V (A3 - Ai) 2 + 46 lC2 772 



It is known that the transition solution (p\ , 0, p* n ) is stable 
near 6c. Therefore the eigenvalues of L satisfy 

0i(T,p) < 0, (3±(T,p) < for (T,p) near cd. 

However, near /ic there is a curve segment /i'c such that 



Figure^ by 

E 1 = {(T,p)eR 2 + I (X 1 ,X 2 ,X 3 )(T,p) = (-,-,+)}, 
£2 = {(T,p)£R 2 + I (A 1 ,A 2 ,A 3 )(T,p) = (-_-)}, 
//^{(T^gM 2 , I (Ai,A 2 ,A 3 )(T,p) = (+,-+)}, 
H 2 = {(T, P ) eR\ I (Ai,A 2j A3)(r,p) = (+,+,+)}, 
Region fbc={(T,p)eR 2 + | (Ai, A 2 , A 2 ,A 3 )(T,p) 
= (+,--,-)}, 

iie^on ed/i' = {(T,p) G | (Ai, A 2) A 2) A 3 )(T,p) 
= (--+,-)}, 

and let the Region 0g' f'ch' be the complement of the sum 
of the above regions. Then the following conclusions hold 
true: 

(1) If (T,p) G Ex, the phase of 3 He is in solid state. 

(2) If (T,p) G E 2 , the phase is in normal liquid state. 

(3) If (T,p) G Region f'bc, the phase is in phase A 
superfluid state. 

(4) If (T,p) G Region cdh' , the phase is in phase C 
superfluid state. 

(5) If (T,p) G Region Og'f'ch', the phase is in phase B 
superfluid state. 

(6) If {T,p) G Hi, there are two regions Ri and R2 
in the state space (pi,p 2 ,p n ) such that, under a 
fluctuation which is described by the initial value 
(x ,yo,z ) in (77p: If(x ,y ,z ) G i?i then the 



phase is in solid state, and if (xq, yo, zq) G R 2 then 
it is in phase A superfluid state. 

(7) If (T,p) G H2, there are two regions Ki and 
K 2 in the state space (pi,p 2 ,p n ) such that if 
(xq, yo, Zq) G Ki then the phase is in solid state, 
and if (xq, yo, zq) G K2 then it is in phase B super- 
fluid state. 



IV. CLASSIFICATION OF SUPERFLUID 
TRANSITIONS 



h'c={(T,p)eR 2 + \ (3i(T,p) = Xi=0}. 



Thus, system ( 18 1 has a t ran sition on h'c, which is called 



the second transition of and h'c is the coexistence 
curve of phases C and B; see Figure [2] 

In summary, with the above analysis and the dynamic 
transition theory, we arrive at the following transition 
theorem: 



Theorem 1 Define a few regions in the PT -plane (see 



In this section, we classify the superfluid transitions of 



(11 1 crossing various coexistence curves. 



First we consider the transitions crossing curve seg- 
ments cd and be in Figure [2J Obviously, we have 

cd={{T,p) eR\ I (Ai,A 2 ,A 3 )(T,p) = (- 0,-)}, 
bc={(T,p) GM 2 + I (Ai, A 2 , A 3 ) (T,p) = (Q, -,-)}. 



Let 

Ai = aici - a 3 \X 3 \, 
Bi = bic 2 - b 3 \X 3 \, 
In fact, it is obvious that A 2 



A 2 = a\C2 - a 2 \X 3 \ 
B 2 = hci - b 2 \X 3 \. 

-B 2 . 



Theorem 2 For the system (11) we have the following 
assertions: 



(1) As (TojPo) S cd, the transition of (11) at (Tq,pq) 
is between the phase B and normal liquid. Further- 
more if B x < 0, then it is Type-I, and if B x > 0, 
then it is Type-II. 

(2) As (Tq,pq) £ be, the transition is between the phase 
A and normal liquid. Moreover, if A x < 0, then it 
is Type-I, and if A x > 0, then it is Type-II. 

Theorem [2] provides conditions for the first transition 
of (111. The following theorem gives sufficient conditions 



for the second transition near fc in Figure 
Obviously, the curve fc is given by 

Jc={(T,p) el} | (A 1 ,A 2 ,A 3 )(r,p) = (+,0,-)}. 

To set up the second transition theorem, we need to 
assume the following conditions. Let e > be small. 
Suppose that 



B-2 > —£, b x a 3 - b 2 ai > -e for (T,p) £ fc, 
and the gap between be and fc is small, i.e., 

\T 2 -T 3 \ = 0(s) V(T 2 ,p)ebc, (T 3 ,p)eTc 
We also assume that 

(ax,cx) = O(e), (b x ,c 2 ,c 3 ,a 3 ,b 3 ) = 0(1), 
a 3 Bi - b 2 A 2 > in fc with A x < 0. 



(21) 

(22) 

(23) 
(24) 



Theorem 3 Under conditions (21) and (22), there ex- 
ists a curve segment fc near fc as shown in Figure 



such that (11) has the second transition from the first 
transition solution (/?f,Q,p*), i.e., (15) has a transition 



from (pi, p 2 , p n ) = in fc, and the transition solutions 
(pi,p 2 ,p n ) satisfy that p 2 > 0. In addition, if (23) and 
\2JX hold true, then this transition is Type-II. 



Physical experiments display that the superfiuid tran- 
sition of liquid 3 He between the normal liquid and su- 
perfiuid phase B is continuous. Hence, it is necessary 
to give the conditions of Type-I transition of (11) at the 
intersecting point C of two curves Ai = and X 2 = 0. 

Theorem 4 Let (To,po) be the point C that X x (Tq,pq) = 
and X 2 (Tq,pq) = 0. Then the transition of (11) at 
(Tq,pq) is Type-I if and only if one of the following two 
conditions hold true: 

(i) A 1 <0,B 1 <0,A 2 = B 2 <0, 

(it) A x < 0,Si < 0,A 2 = B 2 > and A X B X > A 2 B 2 . 

In particular, if the transition is Type-I, then for Ai > 
0, X 2 > near (po, To), there are four types of topological 
structure of the transition on center manifold, which are 
classified as follows: 



(1) This transition is of the structure as shown in Fig- 
ure^jfa), if 

AJBil + X 2 A 2 > and A 2 Ui I + X X B 2 > 0. 



(2) The transition is of the structure as shown in Fig- 
ure^b), if 

Ail-Bil + X 2 A 2 < and X 2 \ A x \ + X X B 2 < 0. 

(3) The transitions is of the structure as shown in Fig- 
ure^c), if 

X X \B X \ + X 2 A 2 < and X 2 \ A x \ + X X B 2 > 0. 

(4) The transition has the structure as shown in Figure 
Wd), if 

X X \B X \ + X 2 A 2 > and A 2 Ui| + Ai-B 2 < 0. 







FIG. 6: Topological structure of Type-I transition near the 
intersection point C = (T ,p ) of Ai = and A 2 = 0. 

Before the proof of these theorems., we have the fol- 
lowing remark. 

Remark 1 Physically, the transition between normal 
liquid and superfiuid for 3 He is generally Type-I, the 
transition between superfiuid phases A and B is Type-II, 
and the region of phase A is narrow. Therefore, under 
the conditions ( 21 1-( 24 ) and 



B x < in cd, 

B x < 0, A x < near point C = f x nf 2 = (T ,p ), 



the above theorems (Theorems 2]|4) provide a precise 
mathematical proof for superfiuid transitions of liquid 
3 He with no applied magnetic field. 
By condition ( 23 1 we see that 



a x -c x = 0{e 2 ), a 3 = 0(l). 

Assertion (1) of Theorem [2] implies that only in a very 
small range of (T,p) near the change point of superfiuid 
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and solid, the transition between normal liquid and su- 
perfluid is II-type, and this range is 

0>A 3 (T,p)>-^ = -O( £ 2 ). 

Moreover, the superfluid density of phase A near the solid 
phase is in the quantitative order e 3 , i.e., 



Pi 



ajc 3 



0(e 3 ). 



Hence, difference between the Type-I and the Type-II 
phase transitions in experiments is very small. 

Note that Theorem [3] are also valid if condition (21 ) is 
replaced by that 



8X2 
dT 



3> 1 for (T,p) near fc. 



Proof of Theorem [2} As (T 0) p ) G bc,\ 2 (T ,p ) < 
and the space of (pi, 0, p n ) is invariant for (111. There- 
fore, the transition equations of (111 at (Tq,po) are re- 
ferred to the following form 



dpi 
dt 

dfn 
dt 



- aip n pi - a 3 p 2 , 

: hPn ~ C1P1 + C 3 pl 



c 4 p 3 n . 



(25) 



The second order approximation of the center manifold 
function p n of (251 satisfies the equation 



A3P77 



C 3Pn 



C\p\. 



Its solution is 

Pn 



C?C3 2 



ClPl 



|A 3 | |A 3 | 3 



Pt + o(pi). 



Putting p n in the first equation of (251 we get the reduced 
equation of (111 on center manifold as follows 

l -Mp\- a -^pl + o{pl). (26) 



dpi 
dt 



I A; 



Assertion (2) follows from (26 1 



Likewise, if (Xo,po) G cd, Ai(To,po) < and the space 
of (0,p 2 ,p n ) is invariant for (111, therefore in the same 



fashion as above we can prove Assertion (1). The proof 
is complete. ■ 

Proof of Theorem [3j We proceed in the following 
two cases. 

Case 1: A\ < in be. In this case, by Theorem [2] 
the transition of ( |ll| in be is Type-I, and the transition 
solution (pi,0, Pn) satisfies that 



Pn 



Cl 



rPi- 



The equations describing the second transition are given 
by ( 15 ) and the eigenvalues in ( 16 1 and ( 17 1 are rewritten 
as 



/?! =A 2 + 61 1/5* I - b 2 p\ 



=A 2 + tt— rB 2 pt, 

l A 3| 



[Ai + A 3 + J(Ai - A3) 2 + 4aiCip*J 



=(by (22) and p\, \\ being small) 
= -Ai + o(|A 1 |), 
Ik </?+■ 
In addition, we know that 



(27) 



(28) 



(29) 



Ai = 0, A 2 < 
X 2 (T-S, P ) > 



on be, 

for (T,p) e fc and 6 > 0. 



Hence, by assumptions (211 and &22h, from (27l-(29l we 



can infer that there exists a curve segment fc near fc 
such that for (T 2 ,p) € be and (To,p) £ T' 3 we have 



<0 if T <T<T 2 , 
k(T,p)\ -0 ifT = T 0) 
>0 ifT<T , 



(30) 



and /3_ < (3 + = - 
[6] the system ( 15 
To determine t. 



Ai -aipl 

-Ci A 3 



-Ai + o(|Ai|) < 0. Hence, by Theorem 

has a transition on T' 3 . 
ie transition type, we consider the cen- 

(31) 



ter manifold function of (|15|, which satisfies 

o(p 2 ) 



Pi 

Pn 



a.2P\P2 
c 2 p 2 



The solution of (31 ) is 

a 2 A 3 + c 2 ai 



Pi 



A3A1 - cxaip\ 
Aic 2 + c\a 2 pX 
A3A1 - Ciaipl 



P*iPi + o(p 2 ), 

P2 + 0{p 2 ). 



From (|16|, (|22| and (|23|) we can obtain 

c 2 



Pi 



A 2 p 2 
a 3 1 A 3 | : 



|A 3 | 



Pi- 



(32) 



Inserting the center manifold function ( 32 1 into the sec- 

3t ■ 
1 



ond equation of (15), we get the reduced equation as 
dp 2 



dt 



= A 2 /9 2 



(u b 2 A 2 2 
A3 03 



By (241, the transition of (15) is Type-II 



Case 2. A\ > in be. In this case, the transition of 
Ilk in 6c is Type-II, and the transition solution in be is 



* a i 1 * 1 

Pi = —\Pn\ 

0-3 



Pn 




a 3 c 2 



- 1 
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The eigenvalue j3\ in (17 1 reads 



Px = A 2 + — (61 a 3 - b 2 a 1 )\p* n \. 

«3 



By (|2l) and ([22] it implies that (|30| holds. Hence 
has a second transition in T' 3 for Ai > in 6c. 



Under the condition ( 23 1 , we have 
a\Ci 



Pn 



03 C3 



Ai 



-«3Pl, 



A3 ~ A3 



1a\C\ 

013 



By A\ > we get that |A 3 | < 0(e 2 ). Thus, the solutions 
of (31 1 can be rewritten as 



Pi 

P2 



aici t c 2 

~~ -PlP2 ^ — P2, 

A1A3 - ciaipj ci 

_03_C2 
O1C1' 



Since ([8} have variational structure, the flows of (34 1 
are of gradient type. Therefore, (pi, p 2 ) = has no ellip- 
tic region for (34). Hence, in the same fashion as used in 
Section 6.3 in |2], one can prove that the region 

S = {( Pl , P2 ) eK 2 | P i >0,p 2 >0} 

is a stable parabolic region. Namely, ( P i,p 2 ) — is an 
asymptotically stable singular point of p4| if and only if 
one of the two conditions (i) and (ii) holds true. Thus, 
we only need to prove Assertions (l)-(4). 

For Type-I transition at point C = (To,po), by con- 
dition (i) and (ii), A\ < and B\ < 0. Hence, as 
Ai > 0, A2 > there are bifurcated solutions of (33 1 
in the pi-axis and /92-axis as 



-P2- 



«i = (pi,0) 



obtain reduced equation on the center manifold as 



Putting p 1 and p n into the second equation of ([15} , we The Jacobian matrices of p3l at z x and z 2 are given by 



df2 
dt 



A2P2 



hc 2 a 3 c 2 b 2 



aici 



Due to ( 23 1 we see that 

b x c 2 a 3 c 2 b 2 



ci 



63 > 0. 



&3 ^- 



a\Ci 



ci 



Therefore, this transition of (15) is Type-II 



It is clear that the second transition solutions 
(f>i,p2,Pn) satisfy that p 2 > 0. Thus, the theorem is 
proved. ■ 

Proof of Theorem [4j At point C = 

(T Q ,p ), \i(T ,p ) = 0, \ 2 (T ,po) = 0. Hence, the center 
manifold function of (ITT]) at (T , po ) reads 



Pn = 



Cl 

I A, 



e-2 



c 3 



x — rP2 + tt^ (ciPi + c 2 p 2 ) 2 - 
A 3 | |A 3 | d 



Putting p rl in the first two equations of (111 we get the 



reduced equations on the center manifold as 

dpi 
dt 



dpi 1 2 

AlPl + 7T— r(^lPl + ^2/Ol/0 2 ) 
l A 3| 



C30l_ 

I A3 1 3 



(cipi + c 2 p 2 ) 2 pi, 



(33) 



<^/°2 1/2 

— =A 2 p 2 + j^j{Bxp 2 + B 2Pl p 2 ) 

C3&1 , . o 
- |A^|3^ ClPl + C2P2 > P2 ' 

To verify the Type-I transition, by the attractor bifur- 
cation theorem [2J, it suffices to consider the following 
equations: 



dpi 
dt 

dp 2 
dt 



A x p\ + A 2 p lP2 - j^-k(ciPi + c 2 p 2 ) 2 p 1 , 
B x p\ + B 2 p x p 2 - tt-^{ciPi + c 2 p 2 ) 2 p 2 . 



(34) 



J(z x ) 



-Ai 




Al + JbV\ X2 







Since ( 33 1 has at most three bifurcated sing 



region S, there are only four types of Type- 
as shown in Figure |6ja)-(d), and each type 
determined by the signs of the eigenvalues 
J{z 2 ). Thus, from the eigenvalues f3(z x ) — 

X X B 2 ) of J{z x ) and 0(z 2 ) = ^(A^l + A 
it is readily to derive Assertions (l)-(4) 
complete. ■ 



;ular points in 
-I transitions, 
is completely 
of J(z x ) and 

2 A 2 ) of J(z 2 ) 
The proof is 



V. LIQUID 3 HE WITH NONZERO APPLIED 
FIELD 

When liquid 3 He is placed in a magnetic field H, the 
superfluid transition is different from that with H = 0. 
Experiments show that as a magnetic field is applied, 
a new superfluid phase A x appears, and the region of 
phase A can be extended to the bottom at p = 0. The 
PT-phase diagram is schematically illustrated by Figure 
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FIG. 7: PT-phase diagram of 3 He in a magnetic field. 

As a magnetic field H is applied, there is a new pairing 
state to appear, in which the spin of pairing atoms is 
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parallel to the magnetic field. This new state corresponds 
to the phase A\, and is expressed by 

V2|$ >= 2a\ TT> ■ 

We introduce the complex valued functions tpo to the 
state | TT>,t/>i to the state | ll>, and ip2 to this state 
I tl> +1 4T>- Let po,p a ,Pb stand for the densities of 
superfluid phases Ai,A,B respectively. Then we have 



Po = iV'ol 2 , 
Pa = Toluol 2 
Pb = T 2 \tp a \ 2 + t^I 2 + 



(r > 0,ri > 0), 

(r 2 ,T 3 >0,T 4 >0), 



and the total density of liquid 3 He in a magnetic field is 
given by 



Pn + Po at state A\ , 
P = \ Pn + Pa at state A, 
p n + pi, at state B. 



Thus, similar to ([lj, for liquid 3 He with H ^ we give 
the Ginzburg-Landau free energy in the following form. 
For simplicity we take the nondimcnsional form: 



G(ip ,ipi,ip2,Pn) =\ I U|V^ | 2 - Aol^ol 2 + aoPnlV'ol 2 + ailV'oHV'il 2 + a 2 MH 2 + ttIV'oI 4 



+ Kl |V^i| 2 - \iHi\ 2 + aiP„|Vi| 2 + a 2 |^i| 2 |V'2| 2 + y l^il 4 

+ ^|VV 2 | 2 ~ A 2 |V 2 | 2 + b lPl \^ 2 \ 2 + + « 3 |Vp„| 2 - A 3 |p n | 



2 _ £3 3 _ C4 4 



dx. (35) 



The equations describing liquid 3 He with H ^ read 

kqAV'o + AqV'o - aopnipo - ai\ipi\ 2 ip ~ c^lV^lVo - a 3 \i>o\ 2 ipo, 



dtpo 

at 

dtpi 
dt 

dtp 2 
dt 

dp n 

dt 
a 



KiAf/'i + Xiipi - aip n ipi - ai\tpo\ 2 ipi - a 2 |V'2| 2 V'i - o-3|'0i| 2 '0i 

K 2 Alp 2 + A 2 ^2 - hp n 1p2 ~ "2|V'o| 2 V'2 ~ «2 | V'l 1 2 ^2 ~ b 3 \l/> 2 \ 2 fa 



\ ^0 i / 1 2 ^1 i 1 2 ^1 i 1 2 2 3 

K 3 Ap„ + A 3 p„ - — l^ol - ylW -7^2 1 + c 2 p„-c 3 p„, 



On 



on dil, 



(36) 



r 



where the coefficients satisfy that for any < i < 3 and mogeneous on ft, the problem (36 1 can be reduced to 
1 < 3 < 3, 



a, > 0, cu > 0, &i,&3,C2,c 3 > 0. 



Equations (36 1 should be the same as fl8| for H = 0. 



Therefore we assume that when H = 0, 



dpi 
dt 
dpi 

k q — kx, Ao = Ai, ao = oi, ai = 0, a 2 = a 2 , a 3 = a 3 . dt 

(37) dp» 

Based on the physical facts, we also assume that dt 



\ po - a p n p - axpxPo - a 2 p 2 po - «3Po, 
Aipi - aipnpi - aip pi - a 2 p 2 p\ - a 3 p 2 , 
A 2 p 2 - bip n p 2 - a 2 poP2 ~ a 2 p 2 pi - 6 3 p 2 , 
■ hp 



a ax bx 23 
y Po - yPi - yP2 + c 2 p„ - c 3 p„, 



\ = \ 1 (T,p) + \(T lP ,H), 
X(T,p,H)>0 if if 7^ 0, 

X(T,p,H) -»0 if if ->■(). 



(38) 



( 39 ) 

where p, = |^| 2 (z = 0,1,2), and A; (i = 1,2,3) are as 
in pi. 



Let Ai,A 2 , and A 3 be that as shown in Figure [3^a)- 
(c) respectively. Then, due to (38 1 the curves Xj(T,p) = 
0(0 < j < 3) in PT-plane are schematically illustrated 
When the magnetic field H and the pressure p are ho- by Figure [8] 
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FIG. 8: Curve agh is Ao = 0, bfh is Ai = 0, cem is A2 = 0, 
and defy is A3 = 0. 

Let the applied magnetic field H ^ such that 

< \{T,p, H) < e, for e > small, (40) 



where A is as in (38) 



Then, under the conditions 
(|21|,(|22|,(|37||,(l38| and (|40j, by using the same fashion 
as in Theorem 2|"a nd|3| we can prove the following tran- 



sition theorem for ( 39 1 



Theorem 5 Assume the conditions (21), (22), [3D, (38) 



and (4-0), then for H 7^ there exist two curve segments 



f'h' near Ai = and e'm' near A2 = in the PT-plane 
as shown in Figure [5| such that the following assertions 
hold true: 



(1) The system (39) has a transition in curve segment 
Ao = with A3 < (i.e., the curve segment gh in 
Figure^, which is Type- 1 for (Xq — 2|A3|«3 < 0, 
and is Type-II for — 2| A3 > 0. The transition 
solution is given by (/9q,0,0, p* n ) with pg > 0, p* n < 
0. " " 

(2) The system has a second transition from 
(Pq,0,0, p* n ) in the curve segment f'h' (i.e., 
Ai — a\p* n — ct\p*Q = 0), and the transition solution 
is as (p' , p[,0, p'n) with p' Q > 0,p[ > and p' n < 0. 

(3) The system has a third transition from 
(p' , Pi, 0, p' n ) in the curve segment e'm' (i.e., 
A2 — b\p' n — «2Po — °2p'i = 0), and the transition 
solution is {p'o,p",p2,Pn) mth Pi > 0(0 < i < 2) 
and p' n < 0. 



Remark 2 The first transition of ( 39 ) in curve segment 
gn corresponds to the phase transition of 3 Hc in a mag- 
netic field between the normal liquid and supcrfluid phase 
Ai, and the second transition in f'h' corresponds to the 
phase transition between superfluid phase A\ and A, and 
the third transition in e'm' corresponds to the phase tran- 
sition between superfluid phases A and B; see Figure [9] 

Remark 3 The transition theorems, Theorems ??-[3]and 
[5] provide theoretical foundation to explain the PT-phase 
diagrams of superfluidity, meanwhile they support these 
models of liquid He which are based on the phenomenol- 
ogy 



VI. PHYSICAL REMARKS 

By carefully examining the classical phase transition 
diagrams and with both mathematical and physical in- 
sights offered by the dynamical transition theory, we de- 
rived two new models for superfluidity of 3 He with or 
without applied field. A crucial component of these 
two models is the introduction of three wave functions 
to represent Anderson-Brinkman-Morel (ABM) and the 
Balian-Werthamer (BM) states. 

Then we have obtained a theoretical PT-phase dia- 
gram of 3 He as shown in Figure [2] based on the models 
and the dynamic phase transition analysis. A few main 
characteristics of the results are as follows. 

First, the analysis shows the existence of a unstable 
region H = Hi U H2, in which the solid state and the 
superfluid states A and B appear randomly depending 
on fluctuations. In particular, in Hi, phase B superfluid 
state and the solid may appear, and in Hi the phase A 
superfluid state and the solid sate may appear. 

Second, theoretical analysis suggests the existence of 
phase C superfluid state, which is characterized by the 
wave function ip2, representing | fj,> +| |T>- However, 
phase C region is very narrow, which may be the reason 
why it is hard to be observed in experiments. 

Third, the curve bed is the first critical curve where 
phase transition between normal fluid and superfluid 
states occur. The curve f'c is the coexistence curve be- 
tween phases A and B. The curve be is the coexistence 
curve between normal fluid state and the phase A super- 
fluid state, the curve cd is the coexistence curve between 
normal fluid state and the phase C superfluid state, the 
curve ch' is the coexistence curve between the phases B 
and C superfluid states. 

Fourth, Theorems [2]|4] imply that near the two triple 
points b and c, there is a possibility of the existence of two 
switch points, where the transition on the corresponding 
coexistence curve switches types at each switch point. 
The existence of such switch points depends on the phys- 
ical parameters. 

In comparison to classical results as shown in Figure [T] 
our results lead to the predictions of the existence of 1) an 
unstable region H, 2) a new phase C in a narrow region, 
and 3) switch points. It is hoped that these predictions 
will be useful for designing better physical experiments 
and lead to better understanding of the physical mecha- 
nism of superfluidity. 




FIG. 9: 
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APPENDIX A: GENERAL PRINCIPLES OF 
PHASE TRANSITION DYNAMICS 

In this appendix, we introduce a new phase dynamic 
transition classification scheme to classify phase transi- 
tions into three categories: Type-I, Type-II and Type- 
Ill, corresponding mathematically continuous, jump, and 
mixed transitions, respectively. 



1. Dynamic transition theory 

In sciences, nonlinear dissipative systems are gener- 
ally governed by differential equations, which can be ex- 
pressed in the following abstract form Let X and X\ be 
two Banach spaces, and X\ C X a compact and dense in- 
clusion. Hereafter we always consider the following non- 
linear evolution equations 



du 



L x u + G(u,X), u(0) = tp, 



(Al) 



where u : [0, oo) — > X is unknown function, and A € I 1 
is the system parameter. 

Assume that L\ : X\ — > X is a parameterized linear 
completely continuous field depending contiguously on 
AeK 1 , which satisfies 

L\ = —A + B\ a sectorial operator, 

A : X\ — > X a linear homeomorphism, (A2) 

B\ : X\ —* X a linear compact operator. 

In this case, we can define the fractional order spaces X a 
for a G R 1 . Then we also assume that G(, A) : X a — > X 
is C r (r > 1) bounded mapping for some < a < 1, 
depending continuously on A £ K 1 , and 



G(u,X) = o(\\u\\ 



x 



VA e 



(A3) 



Hereafter we always assume the conditions ( A2 1 and 



(A3), which represent that the system (All has a dissi 



pative structure. 

A state of the system (Al I at A is usually referred to 
as a compact invariant set H\. In many applications, 
£a is a singular point or a periodic orbit. A state Y,\ of 
( Al ) is stable if Sa is an attractor; otherwise Sa is called 



unstable. 



Definition 1 We say that the system (Al) has a phase 
transition from a state at A = Ao if £a is stable on 
X < Xq (or on X > Xq) and is unstable on X > Xq (or 
on X < Ao )■ The critical parameter Ao is called a critical 
point. In other words, the phase transition corresponds 
to an exchange of stable states. 



Obviously, the attractor bifurcation of (Al I is a type of 
transition. However, bifurcation and transition are two 
different, but related concepts. 



Let {(3j(X) G C | j G N} be the eigenvalues (counting 
multiplicity) of L\, and assume that 

( <0 ifA<A , 

Re&(AW =0 ifA = A , VI < i < m, (A4) 

[ > if A > A , 
Re Pj(X ) < Vj > m + 1. (A5) 

The following theorem is a basic principle of transitions 
from equilibrium states, which provides sufficient condi- 
tions and a basic classification for transitions of nonlinear 
dissipative systems. This theorem is a direct consequence 
of the center manifold theorems and the stable manifold 
theorems; we omit the proof. 



Theorem 6 Let the conditions (A\) and (A5) hold true. 
Then, the system (Al) must have a transition from 
(u, A) = (0,Ao), and there is a neighborhood U C X of 
u = such that the transition is one of the following three 
types: 

(1) Continuous Transition: there exists an open 
and dense set U\ C U such that for any (p G U\, 
the solution u\(t,tp) of \A1^ satisfies 

lim limsup ||u,x(i, f)\\x — 0. 

(2) Jump Transition: for any A < A < A + e with 
some e > 0, there is an open and dense set U\ C U 
such that for any ip G U\, 

limsup || u\ (t,<p) || x >S>0, 

t — >oo 

for some 5 > independent of X. 

(3) Mixed Transition: for any A < A < A + e with 
some e > 0, U can be decomposed into two open sets 
Ui and U2 (Uf not necessarily connected): U = 
U£ + U£ n C/ 2 A = 0, such that 

lim limsup ||u(t, <p)\\x — V<y3 G £7-j\ 

^^-^0 t— »oo 



limsup \\u(t, ip)\\x > S > 

t — >oo 



V(y5 G £/ 2 A 



With this theorem in our disposal, we are in position 
to give a new dynamic classification scheme for dynamic 
phase transitions. 



Definition 2 The phase transitions for \Alty at X = Xq 

is classified using their dynamic properties: continuous, 
jump, and mixed as given in Theorem^ which are called 
Type-I, Type-II and Type-Ill respectively. 

An important aspect of the transition theory is to de- 
termine which of the three types of transitions given by 
Theorem [6] occurs in a specific problem. A corresponding 
dynamic transition theory has been developed recently by 
the authors for this purpose; see [5]. We refer interested 
readers to these references for details of the theory. 
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2. New Ginzburg-Landau models for equilibrium 
phase transitions 

In this section, we recall a new time-dependent 
Ginzburg-Landau theory for modeling equilibrium phase 
transitions in statistical physics. 

Consider a thermal system with a control parame- 
ter A. By the mathematical characterization of gradi- 
ent systems and the le Chatelier principle, for a sys- 
tem with thermodynamic potential TL{u, A), the govern- 
ing equations are essentially determined by the functional 
H(u,X). When the order parameters (u%,--- ,u m ) are 
nonconserved variables, i.e., the integers 



Ui(x,t)dx — cii(t) ^ constant. 



then the time-dependent equations are given by 
on ■ S 

-^=-[3 l —H(u,\) + <t> l (u,Vu,\), (A6) 
at oui 

for 1 < i < m, where Pi > and <&i satisfy 

/ V $~H(u,\)dx = 0. (A7) 
Jn , dUi 



The condition ( A7 1 is required by the Le Chatelier prin- 



ciple. In the concrete problem, the terms $i can be de- 
termined by physical laws and (A7l. We remark here 



that following the le Chatelier principle, one should have 
an inequality constraint. However physical systems often 
obey most simplified rules, as many existing models for 
specific problems are consistent with the equality con- 
straint here. This remark applies to the constraint (A13) 
below as well. 

When the order parameters are the number density 
and the system has no material exchange with the exter- 
nal, then Uj (1 < j < m) are conserved, i.e., 



Uj(x, t)dx = constant. 



(A8) 



This conservation law requires a continuity equation 
nil ■ 

-| = -V-JiKA), (A9) 



where Jj(u, A) is the flux of component uj, satisfying 



jj = -fc 3 v(/ij -^2 fit), 



(A10) 



where /i; is the chemical potential of component ui, 
5 



E 

«t4 



Hi = 



Su 



-H(u,\) - (j}j(u,Vu,X), (All) 



and <pj(u,X) is a function depending on the other com- 
ponents Ui (i j). Thus, from (A9)-(A11 1 we obtain the 
dynamical equations as follows 



dt 



- — H(u, A) — 7 (m, Vm, A) 

OUj 



(A12) 



for 1 < j < m, where (3j > are constants, and <pj satisfy 



5 ii 



-H(u,X)dx = 0. 



(A13) 



When m = 1, i.e., the system is a binary system, 
consisting of two components A and B, then the term 
4>j = 0. The above model covers the classical Cahn- 
Hilliard model. It is worth mentioning that for multi- 
component systems, these <f>j play an important rule in 
deriving good time-dependent models. 

If the order parameters (u\, ■ ■ ■ ,Uk) are coupled to the 
conserved variables • ■ • , u m ), then the dynamical 

equations are 



~dt 
dt 



-PA 



Su 



-7i{u, A) — 4>j (u, Vu, A) 



(A14) 



for 1 < i < k and k + 1 < j < m. Here <&i and <fij satisfy 
(A7l and (A13l, respectively. 



The model ( A14 1 we derive here gives a general form of 



the governing equations to thermodynamic phase transi- 
tions, and will play crucial role in studying the dynamics 
of equilibrium phase transitions in statistic physics. 
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